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ABSTRACT 

The eclipsing and strongly interacting binary star system R Arae (HD149730) is in a 
very active and very short-lived stage of its evolution. R Ara consists of a B9V primary 
and an unknown secondary. We have collected the International Ultraviolet Explorer 
(lUE) archival data on R Ara, with most of the data being studied for the first time. 
There are 117 high resolution lUE spectra taken in 1980, 1982, 1985, 1989, and 1991. 
We provide photometric and spectroscopic evidence for mass transfer and propose a 
geometry for the accretion structure. We use colour scale radial velocity plots to view 
the complicated behavior of the blended absorption features and to distinguish the 
motions of hotter and cooler regions within the system. We observed a primary eclipse 
of R Ara in 2008 and have verified that its period is increasing. A model of the system 
and its evolutionary status is presented. 
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1 INTRODUCTION 

R Ara is an eclipsing binary with a period of just over 4.4 
days and an apparent visual magnitude of 6.7. A secondary 
eclipse has never been clearly seen at any wavelength. Spec- 
troscopically, R Ara is single-lined at best, as the few lines 
that do follow the orbital motion of the primary deviate 
from it for much of the orbit. R Ara's spectrum is plagued 
with blended absorption features and non-orbital motions. 

The Hipparcos data on R Ara indicate a parallax of 
12.44 ± 2.03 mas, which puts it at a distance of roughly 
80 pc. Its proper motion is reported as —4.41 ± 1.75^^|^ in 
right ascension and -13.28 ± 1.49^^ in declination. R Ara 
might be part of the Upper Centaurus-Lupus (UCL) O-B 
association, but on the southern most edge. Since the aver- 
age distance to UCL members is about 140 pc, R Ara would 
be either on the nearest edge of UCL or in the foreground. 
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Also, the proper motions of UCL members are somewhat 
greater than that of R Ara but in the same direction. If a 
member of UCL, R Ara's age could be estimated as roughly 
15 million years. 

R Ara was discovered in 1892 by Alexander W. Roberts 
at Lovendale, at the time it was called "Lacaille 6887", 
or "No. (5949) in Chandler's Catalog". It was reported a s 
"(5949) - Arae ... an Algol-vanahle" by iRobertd l| 18941 ). 
The middle of the 20*'' Century saw more interest in R 
Ara with photome t ric st udies by iHertzsprunel ljl942h and 
iPavne-GaposchkinI (|l945l ) , an d the fir s t (and extremely use- 
ful) spectroscopic analysis bv lSahad3 (|l952|). 

Hertzsprung found the orbital period of R Ara to be 
4"*. 42507 with a primary minimum at JD 2,425,818.028 and 
Payne-Gaposchkin measured a 4*^.42509 period near another 
primary minimum at JD 2,429,433.348. Sahade's spectral 
analysis used 118 plates obtained between 1946 and 1951, 
and he estimated the semi-amplitude of the primary star 
to be Ki ~ 60— and the masses of the components to be 
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Ml = 4M0 and M2 = IAMq. Sahade's analysis determined 
the B9 classification of the primary. 

The Hipparcos data provide support for Sahade's clas- 
sification determination of B9. A star with R Ara's appar- 
ent magnitude, at a distance of 80 pc, would have an abso- 
lute magnitude of -1-2.2, which corresponds to a B9 main- 
sequence star. R Ara's colour of B-V = -0.060 is also indica- 
tive of a B9 star. In addition, the ZAMS model of such a B9 
main sequence star yields a mass of 3.5Mq to 4M0 . 

An original radial velocity (RV) plot bv lSahada l|l952t ) is 
shown in figure [1] These RVs appear to follow orbital motion 
of roughly 60 ^ (Sahade's estimate for Ki) during the last 
quarter of the cycle, although there is no clear orbital motion 
seen in the plot. In Section |4] we provide a Ki = 100 — 
curve on our colour scale RV plots, not as a fit to the data 
but as an aid in comparing and contrasting the observed gas 
motions with estimated orbital motion. 

The first lUE obser va tions were discussed by 
McCluskev fc Kondd (|l983l ). iKondo. et all (|l984l ). and 
Kondo. et al.l (|l985l ) They reported P-Cygni line profiles 



in the 1980 data that exhibit orbital variations, along with 
the presence of C iv and N v indicating high temperature 
regions. Both findings suggest mass transfer processes. Mc- 
Cluskey and Kondo also reported peculiar light variations 
seen in the 1982 lUE data and in the photom etric UV data 
from the Astronomical Netherlands Satellite (|Kondo. et al.l 
Il98ll ). and th e presence of X-rays as seen with the Einstein 
Observatory l|McCluskev fc Konddll983 ). 

The results of these UV analyses spurred more ground- 
based photom etric observations , as reported by groups in 
New Zealand (Nie ld. et aLlll988l : [Forbes, et al.lll98a : iBankj 
I1990I : iNield ,_19_91,) who analysed data taken at Black Birch 
outstation of Carter Observatory. Nield determined an 
ephemeris for the year 1986, for R Ara, of 

Pr.min. = JD 2, 446, 585.1597 -f 4.425132S , 

and she found that a phase shift for primary minimum may 
have occurred. We are very fortunate to have her analysis 
because her ephemeris fits the lUE data better than the 
earlier (1952) ephemeris by Payne-Gaposchkin, and there- 
fore gives us more accurate orbital phase values than were 
used in the analysis of McCluskey and Kondo in the early 
1980s. 

One of Nield's origin al visual light curves is shown in 
figure [1 See iNielH l|l99ll '). pages 246-247, for the parame- 
ters of her model. While her model provides an excellent 
fit to the primary eclipse, the outside-of-eclipse portion of 
the light curve is left essentially un-modeled. Also, Nield's 
model does not allow for Roche lobe overflow (RLOF) and 
therefore doe s not g i ve R Ara the status of an interacting 
system at all. iBanks l|l990l) produced a similar model using 
visual data with an eccentricity-shifted (e = 0.42) secondary 
eclipse and a hot spot to explain the increase in brightness 
before the primary eclipse, but this does not account for the 
other variations, especially those around first quadrature, 
nor does it account for RLOF. It should be noted that such 
high eccentricities are quite unusual for systems with peri- 
ods shorter than 5-7 days because these systems are so close 
that tidal dampening would quickly set the eccentricity to 
zero. While it might be possible for extremely rapid mass 
transfer to drive a non-zero eccentricity, the high eccentric- 
ity reported by Banks is most likely in fact spurious. 



A recent paper by lArias. et al.l (|l998l ) describes some 
optical spectroscopic observations of R Ara taken i n 1991. 
The d escription of the data is similar to that of ISahadd 
(Il952h . characteristic of badly blended absorption and emis- 
sion lines. Arias suggests that the entire system is engulfed 
in a circumbinary envelope. 

This paper focuses on lUE data taken in 1985 and 1989, 
none of which has ever been studied before. These observa- 
tions were well planned as they reveal R Ara in the UV 
during critical parts of its orbit. In section [3] we make sense 
of R Ara's strange photometric variations by modeling it 
with a cool accretion structure surrounding the hot primary 
star. Section |4] tracks the motions of both hot (too hot for 
B9 stars) and cool (too cool for B9 stars) regions. We pro- 
vide the most detailed evidence to date for the existence of a 
transient quasi-disc-like accretion structure delivering mat- 
ter from the less massive secondary star to the more massive 
primary in R Ara. We have detected a hot region on the pri- 
mary that is due to the impact of the mass-transfer stream 
which drives a "hot river" around the primary's equator. 



2 OBSERVATIONS 

The high resolution lUE images are listed in Tables [1] 
through [5l and the orbital phase coverage is also given in 
those tables. The 1989 data are particularly useful because 
there are 36 consecutive spectra spanning one orbital cy- 
cle. Such coverage is not obtainable from the ground, and 
is free from many secular variations which are prominent 
for R Ara. Because these 1989 images are fairly evenly dis- 
tributed through the orbit, they were used to construct UV 
light curves and RV plots. The lUE data from 1985 are 
concentrated near secondary eclipse, another beneficial cir- 
cumstance. The short-wave primary (SWF) camera detected 
wavelengths of 1150 ~ 2000 A while the long- wave primary 
(LWP) and long-wave redundant (LWR) cameras covered 
1850 ~ 3400 A. 

All of the lUE data are available in the NEWSIES 
archive. The lUE data used in this paper were obtained 
from the Multimission Archive at the Space Telescope Sci- 
ence Instute (MAST). STScI is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA 
contract NAG5-26555. Support for MAST for non-HST data 
is provided by the NASA Office of Space Science via grant 
NAG5-7584 and by other grants and contracts. 

We also observed a recent primary eclipse of R Ara in 
2008 in order to check for an ongoing phase shift and appar- 
ent period increase. The observations were made with the 
PR0MPT5 telescope (jReichart. et al.|[2005, 'l . Twelve sets of 
ten, one-second exposures were taken at 30 minute intervals, 
in the V filter. The time of minimum light was calculated 



by the method of lKwee fc Van Wocrde]^ (119561 ). using a pre- 
gram based on one by iGhedinil (|1982l ). The resuhing time 
was HJD 2454541.757 ± 0.002. 



3 UV PHOTOMETRIC ANALYSIS 

In our photometric analysis of the lUE data, continuum 
levels were found for "flat" regions of the spectrum that are 
free of spectral features. These local continua were used to 
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Figure 2. Nield's light curve of visual intensity. This is "figure 5" in the paper bv iNieldl l ll99lf) and was reproduced by permis sion o f 
Kathy Nield. The vertical axis is relative intensity and the horizontal axis is the orbital phase (found using the .Payne- Gaposclikini ill 9451 ) 
ephemeris) in units of degrees. The solid line represents the optimal fit of her model. 
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Figure 3. The lUE light curves of absolute calibrated flux at 1320 A. The 1985 data were taken during two consecutive orbits, rather 
than consectuive images during a single orbit. The filled blue circles are data taken on 25 June 1985 and the open blue circles are data 
taken on 30 June 1985, one cycle later. 



Table 1. High resolution lUE observations from 1980. 



1.2 r 





Observation 


Exposure 


Orbital 


Data ID 


Start Time (UT) 


Time (s) 


Phase 


SWP09026 


5/17/1980 20:09 


1080 


0.074 


SWP09138 


5/27/1980 09:08 


1800 


0.231 


SWP09144 


5/27/1980 22:52 


3000 


0.362 


SWP09151 


5/29/1980 22:20 


3300 


0.809 


LWR07781 


5/17/1980 19:57 


360 


0.071 


LWR07782 


5/17/1980 20:47 


660 


0.079 


LWR07869 


5/27/1980 08:42 


1200 


0.226 


LWR07877 


5/27/1980 22:08 


1800 


0.353 


LWR07888 


5/29/1980 21:42 


2040 


0.802 


LWR07889 


5/29/1980 23:20 


1500 


0.816 



produce our light curves. The local continuum levels near 
1320 A for all of the lUE data are plotted in figure [S] These 
are un-normalised light curves of absolute calibrated flux. 

R Ara exhibits significant photometric variations on dif- 
ferent time scales. The secular variations are quite evident in 
the comparison between 1982 and 1989, for example, as the 
primary minimum reaches a similar value while the system 
as a whole is much brighter in 1982 than in 1989. The 1985 
data show that these variations can also take place from one 
orbit to the very next. These variations confuse light curves 
from ground-based images because they must be folded to- 
gether with data spanning several orbits. Observing from 
space can avoid this by enabling consecutive images from a 
single orbit. 

The 1989 lUE data, being studied for the first time 
in this paper, are extremely useful because they offer good 
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Figure 4. The 1989 lUE light curve of the continuum fluxes 
measured at 1320 A (solid circles) and 2915 A (open circles), and 
normalised to 1.0 at orbital phases 0.809 and 0.817, respectively. 



phase coverage of an entire 4.4 day period with consecutive 
images. These data are free of secular variations. 

Our Ught curves of the 1989 data at 1320 A and 2915 
A are shown in figure |4l In addition to the primary eclipse, 
there are two other "dips" in the light curve. The 1980 and 
1982 data (normalized and plotted together in figure O also 
show both dips. The 1980 observations cover the first dip 
and the second dip was observed in 1982, but together they 
are consistent with the single orbit of 1989. Is one of these 
two dips the second ary eclipse? If so, wh ich one? 

It was noted bv lKondo. et~aLl l|l985h that the second of 
these dips (near phase 0.6) becomes shallower with increas- 
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Table 2. High resolution lUE observations from 1982. 



Table 3. High resolution lUE observations from 1985. 





Observation 


Exposure 


Orbital 


Data ID 


Start Time (UT) 


Time (s) 


Phase 


SWP1714S 


6/08/1982 


07:01 


3300 


0.892 


o vv ir ± ( J.'iy 


6/08/1982 


08:36 


2700 


0.906 


SWP1 71 fi(S 

kJ V V i- -Li J-UU 


6/08/1982 


18:38 


3600 


0.001 


SWP17167 


6/08/1982 


20:39 


4200 


0.021 


SWP1 71 71 


6/09/1982 


07:27 


2700 


0.121 


SWP1 71 72 


6/09/1982 


08:45 


2700 


0.133 


SWP171S^ 

kJ V V i- -Li -LOiJ 


6/10/1982 


18:39 


2400 


0.452 


SWP1 71 Sfi 


6/10/1982 


20:00 


2700 


0.465 


SWP1 71 S7 


6/10/1982 


21:13 


2100 


0.476 


^WPI 71 
O vv ir ± ( J-OO 


6/11/1982 


07:32 


2700 


0.574 


^WPI 71 ^^Q 
o vv ir ± ( J-Oy 


6/11/1982 


08:57 


3300 


0.588 


'^WPI 71 QA 


6/11/1982 


18:43 


2100 


0.678 


SWP 17195 


6/11/1982 


19:55 


1800 


0.689 


SWP17196 


6/11/1982 


20:57 


2400 


0.700 


LWR13439 


6/08/1982 


06:27 


1800 


0.884 


LWR13440 


6/08/1982 


08:00 


1500 


0.898 


LWR13441 


6/08/1982 


09:25 


1500 


0.912 


LWR13446 


6/08/1982 


19:44 


3000 


0.011 


LWR13451 


6/09/1982 


06:57 


1500 


0.115 


LWR13452 


6/09/1982 


08:16 


1500 


0.127 


LWR13469 


6/10/1982 


19:28 


1500 


0.458 


LWR13471 


6/11/1982 


07:01 


1500 


0.567 


LWR13472 


6/11/1982 


08:21 


1800 


0.580 


LWR13477 


6/11/1982 


19:24 


1500 


0.684 


LWR13478 


6/11/1982 


20:29 


1500 


0.694 
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Figure 5. The 1980 and 1982 lUE light curves, normalised as in 
figure |4] The filled circles are relative fluxes at 1320 A, and the 
open circles are at 2915 A. 



ing wavelength. Our new 1989 light curves (figure|4]) confirm 
this and extend the condition to the first dip as well. This 
means that neither of the dips can be the result of the sec- 
ondary eclipse, and both dips are caused by something else; 
they are caused by something cooler eclipsing the primary 
and not the primary eclipsing the secondary. We suspect, of 
course, that a cool accretion structure surrounding the B9 
primary is responsible for these dips in the light curve. What 
is the geometry of this accretion structure that causes the 
variations in the light curve? 

In the attempt to model cool regions (or clouds) located 
in the line of sight to the primary star, we chose to use 
the "spots" feature in the bi nary star modeling prog r am B I- 
NARY MAKER 3 (BM3) by iBradstreet fc SteelmanI (|2002h . 





Observation 


Exposure 


Orbital 


Data ID 


Start Time (UT) 


Time (s) 


Phase 


SWP2fi253 

D VV i- ^ij^ijij 


6/25/1985 


05:32 


2400 


0.394 


RWP2fi255 

kJ VV i- ^ij^ijij 


6/25/1985 


08:10 


2520 


0.419 


RWP2fi256 

kJ VV i- ijiJ^tJU 


6/25/1985 


09:31 


2520 


0.432 


SWP2fi257 

kJ V V -L ^yj^ij 1 


6/25/1985 


10:48 


2580 


0.444 


SWP2fi25S 

kJ VV -L iU.inJCj 


6/25/1985 


12:08 


2700 


0.457 


SWP2fi259 

kJ VV i- ijij^tjy 


6/25/1985 


13:29 


3000 


0.470 


SWP2fi2fiO 

kJ VV i- ^ij^yju 


6/25/1985 


15:37 


2700 


0.490 


SWP2fi2fi1 

kJ V V -L ^yj^yji. 


6/25/1985 


16:57 


2700 


0.502 


SWP2fi2fi2 

kJ VV _L ^ 


6/25/1985 


18:25 


2640 


0.516 


RWP2fi2fi3 

kJ VV i- ijiJiiUiJ 


6/25/1985 


19:43 


2700 


0.528 


RWP2fi31 5 

kJ VV i- ^yjij A-ij 


6/30/1985 


06:06 


3000 


0.530 


SWP2fi31 6 

kJ VV -L ijUiJiU 


6/30/1985 


07:39 


3360 


0.545 


SWP2fi31 7 

kJ V V -L ZjUiJ J_ 1 


6/30/1985 


09:16 


3300 


0.561 


SWP2fi31 8 

kJ VV i- ijUiJiO 


6/30/1985 


10:50 


2700 


0.575 


SWP2fi31 9 

kJ V V i- ijijijiy 


6/30/1985 


12:11 


2400 


0.587 


SWP2fi324 

kJ VV -L ^yjij^'-t 


6/30/1985 


19:07 


2100 


0.652 


SWP2fi325 

kJ VV -L ^yjij^'j 


6/30/1985 


20:15 


1920 


0.662 


RWP2fi643 

kJ VV i- ijiJUrliJ 


9/06/1985 


19:53 


3600 


0.028 


RWP2fi645 

kJ VV i- ijiJUrliJ 


9/06/1985 


22:57 


2400 


0.055 


TjWPnfi2fi7 


6/25/1985 


06:18 


1500 


0.400 


T,wpnfi2fiq 


6/25/1985 


08:58 


1380 


0.425 


Lwpn627n 


6/25/1985 


10:19 


1380 


0.438 


T,WPnfi271 

ij vv -L yjyfij 1 A- 


6/25/1985 


11:37 


1440 


0.450 


ij V v -L yjyj ^ i ^ 


6/25/1985 


12:59 


1500 


0.463 


LWP06273 


6/25/1985 


15:07 


1500 


0.483 


LWP06274 


6/25/1985 


16:28 


1500 


0.496 


LWP06275 


6/25/1985 


17:49 


1440 


0.509 


LWP06276 


6/25/1985 


19:15 


1440 


0.522 


LWP06277 


6/25/1985 


20:32 


1020 


0.534 


LWP06301 


6/30/1985 


05:35 


1500 


0.524 


LWP06302 


6/30/1985 


07:05 


1680 


0.538 


LWP06303 


6/30/1985 


08:41 


1680 


0.553 


LWP06304 


6/30/1985 


10:18 


1500 


0.568 


LWP06305 


6/30/1985 


11:41 


1380 


0.581 


LWP06306 


6/30/1985 


12:57 


900 


0.592 


LWP06311 


6/30/1985 


19:47 


1440 


0.657 


LWP06698 


9/06/1985 


21:00 


1920 


0.036 



We started with the parameters given bv lNieldll|l99lf ). which 
as previously noted were a good fit for the primary eclipse 
but not for the rest of the light curve. Nonetheless, their 
model of R Ara's primary eclipse is an excellent starting 
point for understanding its light curve. 

For our model, we first increased the size of the sec- 
ondary to fill its Roche lobe, and to therefore allow for 
RLOF. This, in turn, led to a slightly greater inclination 
angle. We kept Nield's parameters for gravity darkening, 
fimb darkening, and reflection effects. Table|6]lists the other 
parameters used in our model. Figure |S] shows our synthetic 
fight curves plotted over the 1989 lUE data, for each end of 
the UV spectrum, 1320 A and 2915 A. 

A hot spot, similar to that suggested bv iBank"^ (|l990l ) 
to explain the "hump" preceding primary minumum, has 
effects that are not consistent with the data throughout the 
UV spectrum and it does not account for the other variations 
throughout the orbit. 

Our model's "cool spots" should be interpreted as cool 
clouds in the line of sight to the primary star. Because BM3 
is limited to circular regions located at the photosphere of 
the star, the clouds of our synthetic model have the same re- 
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Table 4. High resolution lUE observations from 1989. 





Observation 


Exposure 


Orbital 


Data ID 


Start Time (UT) 


Time (s) 


Phase 


kj V V -L rjyjij 1 vj 


q/1 n/1 QSQ n7-^^ 


2400 


0.977 


o vv L ouy 1 ^ 


Q /I n/1 QSQ OQ-I 7 


3600 


0.991 


kJ V V -L rjyjcj 1 ij 


Q/1 n/1 Q8Q 1 1 -21 


3000 


0.010 


swpsfitjyfi 

kJ V V -L iJ\J U 1 \J 


Q/1 n/1 QSQ 1 2-^8 


3600 


0.026 


SWPSfiQTY 

kJ V V -L iJ\J U 1 1 


Q/1 n/1 QSQ 1 4-^4 


3600 


0.044 


kJ V V -L rJyJCJ 1 O 


Q/1 n/1 QSQ 1 fi-4^ 


3300 


0.061 


kJ VV -L rJ\JCJ(J\J 


Q/1 1 /1QSQ n^-18 


3000 


0.179 


kJ V V J. iJU C/O 1 


Q/1 1 /1QSQ nfi-^fi 


480 


0.191 


SWPSfiQSS 

kJ VV -L OUCyOCj 


Q/1 1 /1QSQ n7-3^ 


3600 


0.201 


kJ VV -L rJyJCJiJiJ 


Q/1 1 /1QSQ 1 "i-^Q 


3300 


0.280 


kJ V V -L rj\jcju<y 


Q/1 2/1 QSQ n4-nQ 


3600 


0.395 


SWP "57002 

kJ V V J. ij 1 yjyj^ 


Q/1 2/1 QSQ 1 1 -22 


3000 


0.462 


kJ V V J. ij 1 yjyjij 


Q/1 2/1 QSQ 1 f^-W^ 

Zj 1 -Lii/ i-ZjiJZj -LO.iJfJ 


3600 


0.534 


kJ VV -T O ( UUU 


Q /I 9 /I QSQ 90-49 


3300 


0.551 


kJ VV -T O ( UlO 


Q /I /I QSQ 1 n-9fi 


"5000 


u.u / y 


kj VV S: O 1 U 1- i 


Q /I /I QSQ 1 7-"^^ 


3300 


0.747 


kJ VV JT^ O / UZU 


Q/1 A /I QSQ 00-07 


ouuu 




^WP'^7n9'^ 

kJ VV -T O ( U^O 


Q/14./1QSQ 07-0'i 


3600 


0.874 


^WP'^7n9'^ 

kJ VV IT O 1 U^O 


Q /I 4/1 QSQ 1 9-4'i 
yy' -Lt:/ -Lyoy iz.'io 


3000 


0.927 


LWP16314 


Q/10/1QSQ 08-44 


1200 


0.983 


T WPl fi'Hl ^ 


Q/1 n/1 QSQ 1 n-2^ 


2400 


0.001 


LWPIB'^IB 

-LJ V V -1- -LUvJXV^ 


Q/in/1QSQ 12-1Q 


1800 


0.018 


TWPl 7 


Q/1 n/1 QSQ 1 4-n8 


2100 


0.035 


LWP16318 


9/in/1989 16:n3 


1920 


0.053 


LWP16319 


9/in/1989 17:49 


1800 


0.069 


LWP16326 


9/11/1989 n6:16 


1800 


0.187 


LWP16330 


9/11/1989 15:22 


1800 


0.272 


LWP16336 


9/12/1989 n5:22 


1800 


0.404 


LWP16339 


9/12/1989 12:22 


1500 


0.470 


LWP16342 


9/12/1989 18:15 


1920 


0.526 


LWP16343 


9/12/1989 20:03 


1800 


0.542 


LWP16349 


9/13/1989 11:25 


1800 


0.687 


LWP16353 


9/13/1989 18:40 


1800 


0.755 


LWP16355 


9/14/1989 01:15 


1800 


0.817 


LWP16358 


9/14/1989 08:12 


1800 


0.883 


LWP16360 


9/14/1989 13:42 


1800 


0.935 


Table 5. Hig 


1 resolution lUE observations from 1991. 




Observation 


Exposure 


Orbital 


Data ID 


Start Time (UT) 


Time (s) 


Phase 


SWP41934 


6/27/1991 14:08 


1800 


0.053 


SWP41936 


6/27/1991 18:35 


3000 


0.096 


SWP41937 


6/27/1991 20:05 


1920 


0.109 


SWP41955 


6/29/1991 13:30 


1380 


0.498 


SWP41958 


6/29/1991 18:00 


2700 


0.543 


LWP20701 


6/27/1991 14:48 


900 


0.058 


LWP20703 


6/27/1991 19:35 


1320 


0.104 


LWP20713 


6/29/1991 17:20 


900 


0.534 



strictions; however, the model does yield a general geometry 
and temperature for the (accretion) clouds. Figure [7] shows 
the correlation between the BM3 model and the entire ac- 
cretion structure. The cool accretion structure can not be 
a solid, symmetric "disc" or "annulus", as such a structure 
would not cause the phase-dependent variation we see. It 
can not be a constant thickness through time either, or we 
would not see the secular variations that we do. Rather, R 
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Figure 6. The 1989 lUE light curves with our BM3 model. The 
solid circles are data points at 1320 A and the open circles are 
at 2915 A. The solid line is the BM3 model at 1320 A and the 
dashed line is the model at 2915 A. 
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Figure 7. (TOP) An overhead view illustrating the geometry of 
the proposed accretion structure. (BOTTOM) The line-of-sight 
views for the first half of the orbit, showing (left) the cool regions 
from the BM3 model and (right) illustrations of the proposed 
accretion structure. The phase value for each pair of images is 
listed to the left. 
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Table 6. The parameters for our "best-fit" BM3 model for R 
Ara. The relative radius is the size of the star relative to the dis- 
tance between the centres of the stars' masses. The temperature 
factor for a spot (cloud) is its temperature relative to the surface 
temperature of the star. 



Mass Ratio 


q 


0.39 (± 0.05) 


Relative Radii 


r2 


0.192 (± 0.025) 
-1.0 (fill Roche Lobe) 


Ratio of the 
Radii 


li 


0.64 (± 0.05) 


X 1 LlUlldl J_J U. 

minosities 
(at 1320 A) 


T,i 
L2 


n QQR (-1- 1 c; V 1 n~3~i 
0.004 (± 1.5 X 10-3) 


Fractional Lu- 
minosities 
(at 2915 A) 


Li 

L2 


0.874 (± 1.5 X 10-3) 
0.126 (± 1.5 X 10-3) 


Surface 
Temperatures 


Ti 

T2 


12500 (± 2000) K 
7000 (± 1300) K 


Orientation 


Inclination 
Norm. Phase 


78° (± 4) 
0.85 


Spot 1 


colatitude 
longitude 
radius 
temp, factor 


80.2° (± 6) 
139.5° (± 3.5) 

35.7° (± 1) 
0.28 (± 0.05) 


Spot 2 


colatitude 
longitude 
radius 
temp, factor 


94.7° (± 6) 
287.0° (± 3.5) 
40.0° (± 1.5) 
0.30 (± 0.05) 



Ara's accretion structure must be more transient and asym- 
metric. 

The asymmetry of the accretion structure around the 
primary, propose d in our model (t o p of fi gure [T}, is similar 
to that found bv lBovarchuk. et al.l l|2002l ) in their hydrody- 
namic simulations of the interacting system (3 Lyr. Coupling 
this asymmetry with the system's inclination may explain 
the orbital variations in R Ara's Ught curve. The cool clouds 
in the restricted BM3 model would correspond to where the 
accretion structure is closest to the primary and is therefore 
in the line of sight to it, while the farther regions of the cloud 
descend below our line of sight to the star revealing more of 
its hot equator. The alternative explanation that would not 
rely on the inclination and asymmetry would be that there 
are higher density concentrations in the accretion structure 
at the location of the "spots" in our BM3 model. 

The transient nature of the accretion structure, pro- 
posed by our model to account for secular variations, is con- 
sistent with the findings of [Richards fc Ratfil l| 19981 ). They 
found that Algol-type binaries with periods shorter than R 
Ara's, such as f3 Per (P = 2.9 days), exhibit mass-transfer 
streams that strike the accreting star more directly and do 
not form stable accretion discs. At the same time, they find 
that the mass transfer in Algols with periods longer than 
R Ara's, like TT Hya (P = 6.9 days), does not create a 
direct impact and therefore quickly produces a stable accre- 
tion disc. An interacting system like R Ara, with a 4.4 day 
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Figure 8. This is our synthetic BM3 light curve under two further 
conditions. The solid line is our model with cool clouds eclipsing 
the primary, as it would look in the visual (5500 A). The dotted 
line is our model in the ultraviolet (1320 A), but with no cool 
clouds. For the dotted line, we disabled only the "spots" feature 
in BM3, leaving all other parameters the same. 



period, might find itself varying between a "(3 Per-like" state 
and a "TT Hya-hke" state. 

We propose that the accretion rate onto the primary in 
R Ara is variable. During times of rapid accretion, the trans- 
fer stream strikes the primary directly and produces a thin- 
ner accretion structure. This should correspond to brighter 
epochs, like 1982. As the accretion structure builds up, dur- 
ing times of slower accretion onto the primary, it will "push" 
the mass-transfer stream away from the star causing a less 
direct impact and a thicker accretion structure. This latter 
case will cause the system to appear dimmer overall, like in 
1989. As the accretion rate increases, the stream may return 
to the more direct impact state and the system will become 
brighter again. 

For completeness, we have included synthetic light 
curves of our BM3 model under two further conditions, in 
figureO The dotted line is our model, at 1320 A, but without 
any cool clouds (or "spots" in BM3) eclipsing the primary. 
This is essentially what Nield's model would look like in the 
UV. Note that the secondary eclipse should almost vanish 
at 1320 A. The lUE data tell us that secondary "dips" in 
the light curve are in fact deeper in UV. 

The solid line in figure [8] is our BM3 model, cool clouds 
included, as viewed at 5500 A, in the visual. Indeed, our 
model can be extended into the visual region to explain the 
"hump" before primary minumum, the wide and apparently 
shifted "secondary eclipse", and the other variations near 
first quadrature. 

Our observation of the 16 March 2008 primary eclipse 
of R Ara confirms the ongoing period increase detected by 
Nield. In section \S\ we estimate R Ara's rates of period 
change and mass transfer. 

We have proposed a model for R Ara to explain its 
strange variations in brightness. In the following section we 
put our photometric model to the test through our analysis 
of R Ara's spectral features. Are the motions of the cooler 
ions consistent with our proposed accretion structure and 
are the motions of the hotter ions consistent with our pro- 
posed "hot river" driven around the primary's equator by 
the impact of the mass-transfer stream? 
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4 UV SPECTROSCOPIC ANALYSIS 

In section |3] we proposed a model for R Ara consisting of an 
undetected F to K type secondary undergoing RLOF and 
donating mass to the more massive B9 primary. The mass- 
transfer stream strikes the primary, either directly or at more 
of a glancing angle, which drives a "hot river" around the 
primary's equator. The imact site and equator are where we 
would expect to find the hotter ions, like N v and C IV, 
which require much hotter (> 30,000 K) plasmas than are 
typically found in late B-type stars. The ions that are too 
cool for a B9, such as Ni ll and Fe ll, should be found in the 
accretion structure. 

The first spectroscopic test of our model is to see if 
the hot ions are visible only when the primary's equator is 
not eclipsed by the accretion structure. Figure |9] shows the 
profiles of the N v (1242.80 A) and C iv (1548.20 A) lines 
at four important phases during the 1989 observations. The 
first spectrum (taken at phase 0.201) shows very weak N 

V and C iv absorption, if any at all, and the correspond- 
ing geometry of the model infers that a cool cloud hides 
the equator of the primary from our view. The next spec- 
trum, phase 0.395, indicates strengthened N v and C iv, 
labelled on the graphs with arrows, corresponding to where 
the accretion cloud dips below our line of sight revealing the 
primary's hot equator. Similarly, the hot equator is covered 
by the cooler accretion cloud near phase 0.534, in the third 
spectrum of figure[9l and the hot ions disappear. The fourth 
spectrum (phase 0.809) shows very strong N v and C iv ab- 
sorption, as we are looking at the region of impact between 
the gas stream and the star and/or accretion structure. 

Another test is to compare the strengths of the hot and 
cool absorption lines throughout the orbit. Although the ab- 
sorption features are badly blended, we did compute equiv- 
alent widths (EWs) by direct integration of the normalised 
data. Figure [TO] displays EW versus orbital phase for the N 

V and Ni ll absorption profiles, from the 1989 lUE data. 
Again, the 1989 images cover a full cycle and are ideal for 
studying orbital variations. Just after primary minimum, 
around phase 0.05 to 0.1, the strength of the N v line de- 
creases while that of the Ni ll reaches a maximum. This 
is the phase when the accretion structure is most optically 
thick in our line of sight to the primary. Later, at around 
phase 0.4, the N v line strengths reach a maximum while the 
Ni II is quite weak. This, according to our model, is when 
the hot equator of the primary is exposed to us and the cool 
accretion structure is mostly below our line of sight. 

There is also some spectroscopic evidence for our expla- 
nation of R Ara's secular variations. During times of more 
direct mass transfer, the impact of the stream may cause a 
splashing of material over the primary. This would behave 
like an asymmetric "wind", seen when the splashing mate- 
rial is projected onto the hot star. Since a stellar wind is 
detected as a P-Cygni line profile, any asymmetric P-Cygni 
lines may be caused by this splashing. The P-Cygni features 
of 1980, noted by McCluskey and Kondo, appear at phases 
0.07 and 0.23 but not at any other phases. During these 
phases we are looking at the system from behind the sec- 
ondary, with a view from above the impact region. There 
are no P-Cygni features in the 1989 data even at similar 
phases. This is a secular spectroscopic variation supporting 
our explanation of the secular photometric variations. 
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Figure 9. N V at 1242.80 A and C IV at 1548.20 A & 1550.77 A, 
shown along with the modeled line-of-sight geometry for phases 
(A) 0.201, (B) 0.395, (C) 0.534, and (D) 0.809, from the 1989 lUE 
data. The arrows indicate the absorption features in the spectra. 
We only see N v and C IV absorption when the cool cloud is not 
eclipsing the primary's hot equator. 
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Figure 10. Equivalent widths of the N v (1242.80 A) and Ni 11 
(1454.84 A) absorption lines, from the 1989 WE data. The solid 
circles are for N v and the open circles are for Ni II. 



An RV analysis provides the most rigorous test of our 
photometric model by answering the question asked at the 
end of section [S] Due to the complicated blending of the 
absorption lines in R Ara's spectrum, we employed colour 
scale RV plots to visualize the dynamic behavior of the ions 
within the system. Such plots are displayed in figure [TTI for 
the hotter ions, and figure [121 for the cooler ions. The 1989 
data was used due to the good phase coverage over a single 
orbit. The arrows to the right of the colour scale RV plots 
indicate actual data points with a linear interpolation of the 
data between them. The dotted white line in these plots is 
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not a fit to the data, but it is the estimated motion of the 
centre of mass of the primary star {k — 100 km s~^) and 
it is plotted to aid in distinguishing orbital from non-orbital 
motion. 

The "hottest" ion, N v, clearly follows the orbital mo- 
tion of the primary star during the last quarter of the orbit, 
with the obvious deviation of the large shortward shifted 
"flare" near phase 0.4. Si iv also shows a clear orbital com- 
ponent during the last quarter of the orbit, and it exhibits 
the flare before the secondary eclipse. This flare will be dis- 
cussed later with the 1985 data, as they are concentrated 
around this portion of the orbit and provide better time 
resolution allowing us to "zoom-in" to that feature. 

The "cooler" ions exhibit a sort of "S-shape" on the 
colour scale RV plots, which is characteristic of an accre- 
tion disc. They show high receding velocities before pri- 
mary eclipse followed by high approaching velocities after 
the eclipse. Also, they lack the mid-orbit flare that the hot- 
ter ions present. Si ll and Fe ll are cool enough to be found 
in the interstellar medium (ISM), and these colour scale RV 
plots clearly show their interstellar components as vertical 
lines. 

The differences between the motions of the hot and cool 
regions within R Ara yield strong support for our model. We 
expect the cooler ions to follow the motion of a disc or cloud 
moving around the primary. We also expect the hottest ions 
to be present at the primary's equator, as the equitorial hot 
river is driven by the impact of the mass-transfer stream. 

We will now zoom in to the flare seen in the RV plots 
for the hotter ions. Specifically, we will consider the Si iv 
line profiles of 1985, which actually display a clear splitting 
of the line into two blended components. Si iv always gives 
us the strongest, sharpest lines, probably due to its high 
abundance. In figure 031 we diplay the colour scale RV plot 
for Si IV in 1985, focused on phases 0.4 to 0.6. Here we 
see the fiare in much more detail. The observations labelled 
a through / in the figure are modeled with two Gaussian- 
shaped absorption lines blended together. The figure also 
shows the geometry of the system at the time of each image, 
including the "cool clouds" from the BM3 model. 

According to our model, we are looking at both the 
accretion structure (which is descending below our line of 
sight to the centre of the primary) and the primary's equator 
(in the centre of our line of sight) during phases a to / of 
figure [131 Figure [141 tracks the RVs of the two components 
through the 1985 observations. The blue points in the plot 
represent the Gaussian component modeled in blue in figure 
1131 and the red points are for the red Gaussian. The red 
component shows RVs that vary as we would expect the 
accretion disc to vary and is likely arising from (the warmer 
parts of) the accretion structure, which appears to oscillate 
between receding and approaching. The blue Gaussian in 
the figure moves steadily to higher and higher approaching 
velocities in the progression from a to /, and it is most likely 
arising from the star's hot equator as we are looking more 
and more "down the barrel" of the gas stream at each of 
those phases. In order the have a hot region of the primary's 
photosphere reach velocities of 800 km s~^, there must be 
a mass transfer stream driving it. 

Our final spectroscopic test of the photometric model 
is to compare the absorption line strengths and RVs in 1982 
with those in 1989. Figure [15] displays colour scale RV plots 
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Figure 14. The radial velocities of the two components of the 
blended Si IV Une in figure [T51 The blue points are the RVs of the 
blue Gaussian and the red points are for the red component. 



for N V and Si iv, the hotter ions, as well as for Si ll and Fe 
II, the cooler ones, in 1982. Since the 1982 observations are 
fewer in number and not as frequent as in 1989, these plots 
are not interpolated, rather they consist of individual strips 
for each image. The first thing we notice is that the cooler 
ions are much weaker in 1982 than in 1989. In fact, mostly all 
of the Si II and Fe ll absorption is by the ISM in 1982. This is 
further proof for our explanation of the photometric secular 
differences between 1982 and 1989. Recall that the system 
is brighter in 1982 than in 1989, resulting from more direct 
mass transfer and a thinner accretion structure. In addition 
to seeing less of the cool accretion structure, we would expect 
to see more of the primary's hot equator in 1982. Indeed, the 
absorption by the hotter ions is stronger in 1982, and those 
ions are present for more of the orbit. We can even see a 
component of the N v line following the estimated orbital 
motion of the primary during every observed phase in 1982. 
In 1982, we are certainly seeing more of the hot equator 
of the primary's photosphere and less of the cool accretion 
structure than in 1989. 



5 DISCUSSION / EVOLUTIONARY STATUS 

Interacting binaries consisting of a cool (F-K) secondary do- 
nating mass to a hotter (B-A) and less evolved primary by 
RLOF are considered to be Algols. But, as is pointed out by 
[Albright fc Richardsl l| 19961 ). Algols are in a slow and steady 
phase of mass transfer (M ~ 10~^^ to lO"'^ Mq yr~^), which 
is why they are used as laboratories for studying accretion 
processes in binary stars. Is R Ara an Algol? If not, in what 
stage of its evolution does it currently reside? 

We have explained R Ara's photometric variations and 
its complicated spectrum as the result of unsteady accretion 
via a variable accretion structure. We can estimate the mass 
transfer rate using the new 2008 PROMPT observations. We 
will use Nield's ephemeris to calculate the period change rate 
and mass transfer rate between 1986 and 2008. The calcu- 
lated time of primary minimum is C = HJD 2,454,541.547 
and the observed time is O = HJD 2,454,541,757. The num- 
ber of orbits between Nield's ephemeris and the 2008 eclipse 
is n = 1798 and the period from the Nield ephemeris is P = 
4.425132 days. The rate of period increase between the two 
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Figure 11. These are colour scale RV plots for N V, C IV, and Si IV in 1989. These images show the motions of the hotter regions of R 
Ara. The arrows indicate actual data and the white dotted curve represents the estimated 100 orbital motion of the centre of mass 
of the primary star. 
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Figure 12. These are colour scale RV plots for Si II, Ni II, and Fe II in 1989. These images show the motions of the cooler regions of R 
Ara. The arrows indicate actual data and the white dotted curve represents the estimated 100 orbital motion of the centre of mass 
of the primary star. 
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Figure 13. The "Splitting" of Si IV absorption line profiles in 1985. The colour scale RV plot on the left is similar to those in figures 
Illl and ll2l but focused between phases 0.4 and 0.6. The time resolution of the 1985 observations is very good, allowing us to track the 
two components as blended Gaussians. The component corresponding to Si IV in the accretion structure is diplayed in red and the Si 
IV located at the primary's hot equator is in blue. The two components blend to give the green profile which matches the data. The 
line-of-sight geometry from the BM3 model is given for each image, and the phase value is printed above each line profile. 



epochs is then: 

p ^ 2(0 - C) 
n?P 

= 2.94 X 10" 



(± 2.80 X 10-^°) '^"-y' 



day 



(1) 
(2) 



From this we find P/P = 2.43 X 10"® yr'^ , which corre- 
sponds to a timescale of roughly 4 x lO'' years for this phase 
of R Ara's evolution; a very short-lived stage. 

If we assume the masses of the stars to be Mi — 4Mq 
and M2 = I.4M0 (these are Sahade's original estimates, 
but they are also consistent with the mass ratios used both 
in Nield's photometric model and our photometric model in 
this paper), and if we assume conse r vative mass transfer, 
we can use the equation from iKwed I 19581 ) to get a mass 
transfer rate of: 



M 



PM1M2 
3P(Mi - M2) 

1.74 X 10"*^ (± 1.66 X 10"*) ^ 

yr 



(3) 
(4) 



which is much more rapid than typical Algols. In fact, at 
this rate of mass transfer the secondary would lose half of 
its current mass in just a few hundred thousand years. 

We believe R Ara is in the pre- Algol stage of its evo- 
lution, just past the reversal of mass ratio, which would 
be characterised by rapid mass transfer and an unstable, 
transient accretion structure. We project that R Ara's mass 



transfer rate will eventually slow down as it enters a more 
stable state. At that time, R Ara may become a classic Algol 
type binary star. 

Many frequent times of primary minima are desired in 
order to more accurately track a change in the orbital pe- 
riod of R Ara and to ultimately conduct a complete 0-C 
analysis. Spectroscopic data for use in constructing Doppler 
tomograms would also be of great benefit to visualize R 
Ara's mass flow and accretion structure. 
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Figure 15. These are colour scale RV plots for N V, Si IV, Si II, and Fe II in 1982. These data are not dense enough to be interpolated, 
so they apear as individual strips for each image. The white dotted curve represents the estimated 100 orbital motion of the centre 
of mass of the primary star. 
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